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Abstract. We solve the conformal bootstrap equations of the four fermi 
model or 0(N) Gross Neveu model to deduce the fermion anomalous di- 
mension of the theory at 0(l/iV 3 ) in arbitrary dimensions. 
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1 Introduction. 



The observation that in the neighbourhood of a phase transition in a (renor- 
malizable) quantum field theory there is a conformal or scaling symmetry 
has proved to be a useful technique in solving models within some expansion 
scheme, [l], ^, |3|. For instance, one can determine the critical exponents of 
all the Green's functions and these govern the quantum properties of the 
field theory. Also, Polyakov demonstrated that the form of the propaga- 
tors and the 3-point vertices are fixed up to constants by demanding that 
the theory is symmetric under the conformal group, Jl. Consequently, with 
this knowledge various authors formulated a method known as the conformal 
bootstrap to solve models such as (p 3 theory in d-dimensions order by order in 
some perturbative coupling constant, 0, ||, |7|. The cf-dimensional exponents 
could then be compared with the e-expansion of similar exponents in other 
models to ascertain whether they both lay in the same universality class, 
|0]. This is important since for physical systems it is useful to understand 
which class of theories underly and describe the experimental observations 
and subsequently allows one to make further predictions from the theoretical 
model. Following the earlier work of || the ideas surrounding the confor- 
mal bootstrap programme were extended to the O(N) bosonic a model in 
H where 1/N replaced the perturbative coupling constant as the expansion 
parameter with iV large. The critical exponent r] of the fundamental field 
was deduced at 0(l/N 3 ) by solving a set of bootstrap equations which were 
derived by using Parisi's method, [H], since the a model also involves a 3- 
point vertex. One considers the Dyson equations of the 2-point function and 
3- vertex with dressed propagators and vertices. Indeed the 0(l/iV 3 ) calcu- 
lation of H built on earlier 0(1/N 2 ) work of [0, [hJ where the model was 
solved completely at this order by similar conformal techniques. In particu- 
lar, the skeleton Dyson equations for the fields were solved by considering, 
by contrast, the situation of dressed propagators but undressed vertices. It 
was the possibility of solving theories in the large N method this way which 
suggested that the full conformal bootstrap could be exploited to probe the 
model further. 

In this paper, we consider the O(N) Gross Neveu model, [11], and de- 
velop the analogous formalism to deduce the anomalous dimension, rj, of the 
fundamental fermion at 0(l/iV 3 ). The viability of such a calculation was 
also suggested from the 0(1/N 2 ) work of [12] which used the same tech- 
niques of [n| in solving the model at the d-dimensional critical point of 
the theory defined as the non-trivial zero of the /3-function where one has a 
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conformal symmetry. Although the 0(N) Gross Neveu model describes the 
dynamics of the four fermi interaction the fermionic nature of the fields does 
not prevent one from using conformal methods and the four point interac- 
tion can be rewritten as a 3-point vertex plus a bosonic auxiliary field which 
means it has the same underlying </> 3 type structure which is fundamental 
in this area. Indeed one of the motivations for studying such models lies in 
its relation to gauge theories. For instance, if we wish to probe four dimen- 
sional theories with matter fields to high orders in large TV" to discover the 
structure of the quantum theory we must fully comprehend the way in which 
fermions have to be dealt with in the formalism. The 0(N) Gross Neveu 
model provides an excellent laboratory for testing such ideas since it has a 
3-point interaction, in its auxiliary field formulation, which from the critical 
point of view mimics the QED interaction. Essentially one is not obstructed 
by the tedious complications which occur with the inclusion of a 7-matrix 
at the vertex of integration. Moreover, an analysis of the four fermi model 
is necessary and interesting in its own right. For example, there has been 
a resurgence of study of the interaction through its relation to providing an 
alternative mechanism of mass generation in the standard model. In (l^, [l4|] 
it was discussed how the effect of the Higgs boson could be replaced by a 
composite bosonic field built out of the binding of two fermions. Further, 
as our calculations will be in arbitrary dimensions we will be able to deduce 
new and useful information on the structure of the three dimensional model 
explicitly. This is currently of interest since four fermi interactions are be- 
lieved to contribute to models describing high T c superconductivity. Thus 
it is important to have the basic formalism in place in terms of the values 
of the underlying three and four loop Feynman diagrams to ensure that one 
can eventually extend the present work to gain predictions for the realistic 
models which will involve coupling to a U(l) gauge theory. 

The aim therefore of this paper is to provide an extensive analysis of the 
four fermi theory at 0(1/N 3 ) in arbitrary dimensions. In p5[ the conformal 
bootstrap equations which we will solve were derived and checked to ensure 
that one could correctly recover the known 0(1/N 2 ) result of Its ex- 
tension to the next order involves the tedious evaluation of several Feynman 
graphs. Unlike the O(N) bosonic a model, []|], there are fewer graphs to 
consider since graphs with fermion loops with an odd number of fermions 
are trivially zero. However, the structure of fermionic massless Feynman 
diagrams means that they can in several instances be much more difficult 
to compute than their bosonic counterparts. 

The paper is organised as follows. In section 2, we present the basic 
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background to the O(N) Gross Neveu model and introduce the conformal 
structure of the fields of the model which are fundamental to the conformal 
bootstrap programme. The master equation whose solution will yield the 
critical exponent r/ at 0(l/iV 3 ) is derived from the three conformal bootstrap 
equations of the model in section 3 where we also introduce Polyakov's 
conformal triangle approach to formulating the equations. In section 4, we 
discuss in detail the techniques required for the evaluation of the massless 
graphs contributing to the vertex function which contains contributions to 
r\ at 0(l/iV 3 ), in a particular limit of the two regularizing parameters which 
are present. The main result of the analysis is presented in section 5 together 
with several concluding remarks. 



2 Review of the model. 

In this section we recall previous conformal approaches to solving the 0(N) 
Gross Neveu model in the large TV approximation in arbitrary dimensions. 
First, we define the lagrangian we use as, [p"l|], 

L = ~W + ~ Y (2-1) 

where ip l is the fermion field, 1 < i < iV, and 1/N will be our expansion 
parameter. The bosonic field a is auxiliary and eliminating it through its 
equation of motion yields the four fermi interaction explicitly. It is more 
appropriate though to use the formulation (2.1) since the integration rules 
to compute massless Feynman integrals are easier to apply to the situation 
with 3-vertices. In the solution of the model in the usual large N expan- 
sion [11] the a" field becomes dynamical in the true vacuum of the quantum 



theory, though classically it is non-propagating. The coupling constant, g, 
is dimensionless in two dimensions, where the model is asymptotically free, 
[[nfl . Indeed the three loop /^-function of (2.1) has been deduced pertur- 
batively in MS using dimensional regularization and in d-dimensions it is, 



(3(g) = (d-2)g-(N-2)g 2 + (N-2)g 3 + 1 I (N-2)(N-7)g 4 (2.2) 



where the original two loop calculation was carried out in |18[. The result 
(2.2) is what one obtains before setting d = 2 in minimal schemes and 
there is therefore no (i-dependence in the higher order coefficients of the 
coupling constant. Indeed (2.2) is the starting point for analysing (2.1) in 
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the conformal approach in (i-dimensions. For d > 2 one observes that there 
is a non-trivial zero of the /3-function at 



9c ~ (2.3) 

to one loop where d = 2 + e. This corresponds to a phase transition in the 
theory which has been widely discussed and analysed in three dimensions in 
recent years, [19-23]. 

In statistical mechanics it is well known that in the neighbourhood of a 
phase transition physical systems exhibit special properties. For instance, 
various (measurable) quantities display a scaling behaviour in that they de- 
pend only on a characteristic length scale such as a correlation length raised 
to a certain power known as the critical exponent or index. This exponent 
totally characterizes the properties of the critical system. In a completely 
analogous way, in the neighbourhood of a phase transition of a (renormaliz- 
able) quantum field theory one observes that there is a conformal symmetry 
present. In other words Green's functions are massless and obey a simple 
power law behaviour. Further, the exponent of certain Green's functions 
such as the propagator can be related to the appropriate critical renor- 
malization group function through a simple analysis of the renormalization 
group equation. (See, for example, [p4|| .) Indeed the exponent will be a 
function of the spacetime dimension d and any other internal parameters 
of the underlying field theory, which for (2.1) will be N. (At criticality the 
coupling constant is not independent but a function of d and N.) 

Therefore with these simple observations it was possible to solve (2.1) 
in a conformal approach which allowed one to go beyond the leading order 
since the masslessness of the problem simplifies the intractable integrals, 



[12 1, which would otherwise occur in the conventional large N approach 
where the fields are massive. We now recall the essential features which are 
relevant for the 0(l/iV 3 ) calculation we discuss here. First, since we will 
now work in the neighbourhood of g c given by (2.3) in d > 2 dimensions we 
write down the most general form of the propagators of the theory consistent 



with Lorentz and conformal symmetry as, 1121 



Ai B 



as x — > in coordinate space. We have chosen a propagating a field since 
we are solving the theory in the true vacuum. Whilst it is more convenient 
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to compute in x-space one can easily deduce the momentum space forms of 
(2.4) through the Fourier transform 



(x 2 Y 



„ikx 



a(a) 



(2.5) 



valid for all a where we set the spacetime dimension d to be d = 2/j for later 
convenience. The quantity a(a) is defined by a(a) = T(fj, — a)/T(a). In (2.4) 
the quantities A and B are the amplitudes of the respective fields and are x- 
independent, whilst we have defined the dimensions of the fields to be a and 
p. Their canonical dimensions can be determined by a dimensional analysis 
of the action which is dimensionless. However, quantum mechanically these 
dimensions do not remain as their engineering values due to quantum fluc- 
tuations such as radiative corrections. To allow for this scenario we define, 

a = ii+\<q , j3 = I-77-2A (2.6) 

where 77 is the anomalous dimension of the fermion and from an examina- 
tion of the 3-vertex of (2.1), 2A is its anomalous dimension. Both 77 and 
A are 0(1/N) and depend on d and N. In a critical point analysis of the 
renormalization group function, 77 corresponds to the wave function renor- 
malization and 77 + 2A to the mass anomalous dimension. They have both 
been determined at 0(1/N 2 ) in arbitrary dimensions as, [12, 25|, 



m 

Ai 



2r(2/j - !)(// - l) 2 

r(2-/z)r( M + i)r20u) 

2(M-1) 



(2.7) 
(2.8) 



m 



A, 



2(/x-l) 2 



A* 



2(M-1) 2 



+ 3/i + 4(/i - 1) + 2(fi - l)(2/i - (//) 



(2.9) 



3/i(/i-i)e + (2/i-i)* 



(2/x-l)(/i 2 -/ J -l) 



(m-1) 

(2.10) 

where 7/ = £,=i A = ^=1 &i/N\ = ^(2/j-l) - ^(1) + ^(2-/i) 

— ^(/i), 0(/j) = ^'(/-t) ~~ ^'(1) an d i/j(x) is the logarithmic derivative of the 



r-function. The result (2.7) has been given originally in [26] whilst (2.8) 
has been determined independently in ||21|| . It is the aim of this paper to 
deduce 773 in arbitrary dimensions. It is worth noting that the 0(l/iV 2 ) 
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correction to the critical exponent A which relates to the /3-function of (2.1) 
has recently been determined too as, O, [27[|, 



Ai = - (2n-l) m (2.11) 



(m-i) 



2 (2/x-3)/i 
(/i - Wm ~ fa - 2) 



2\ 



t ^ 1 1 n 2 3 1 3 
+ * 7 rrr + — r - 2/T - 

V(/i-2) 2 2(/i-2) 2 2/i 

3uG / 6 \ 2 5 

17 11 1 

+ 



V 2(/i-l) (/i-1) 2 ' 4(/i-2) 2(/x-2) 2 J (2 ' 12) 

where $ = ip'(2[i - 1) - ^'(2 — /i) — V'O-O + i/>'(l) which means (2.1) 
has been solved completely at 0(1/N 2 ). It is important to recognise that 
these analytic expressions have been deduced merely from the ansatze (2.4) 
and the solution of the critical Schwinger Dyson equations, which therefore 
means that the method of |], [lOj, [l^] is an extremely powerful one. 

With (2.4) it is possible to deduce the scaling behaviour of the respective 
2-point functions which we will require here. Their coordinate space forms 
are deduced by first mapping (2.4) to momentum space via (2.5), inverting 
the propagator before applying the inverse map to coordinate space, M. 



Thus, we have, [12| 



where the functions have an analogous scaling structure but the amplitudes 
are non-trivially related. The functions r(a — 1) and p(0) are defined as 

m - . Ha) = 1*!SL (2.14) 

It is the quantities (2.13) and (2.14) which will be necessary to solve the 
conformal bootstrap equations. 

3 Conformal bootstrap equations. 

In this section, we recall the basic bootstrap equations for the model which 
will then be solved at 0(1/N 2 ) for r\ and which were given in [15| for (2.1). 
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The development of a bootstrap programme to solving field theories dates 
back to the early work of [||, || . Basically one simplifies the set of Feynman 
diagrams contributing to a Green's function by first using dressed propaga- 
tors. Further, each vertex of the Feynman graph is replaced by a conformal 
triangle which was discussed originally in j|] and later in . To understand 
this concept we first have to introduce the technique known as uniqueness 
which was first discussed in |2^] and is a fundamental tool for evaluating 
massless Feynman diagrams. The rule for a purely bosonic vertex was given 
in [Eq, PUT] whilst the development to the style of vertex which arises in the 



model (2.1) was given in [12] and is illustrated in fig. 1. When the arbitrary 
exponents are constrained to be their uniqueness value Ya=i A = 2/i + 1 
then one can compute the integral over the vertex of integration of the left 
side of fig. 1 to obtain the product of three propagators where fh, fa) 
= tt^ nf=i a (Pi)- The origin of this uniqueness condition can be understood 
in several ways. First, if one carries out the explicit integral using Feynman 
parameters over the vertex of integration one finds that for general the 
calculation cannot be completed due to the appearance of a hypergeometric 
function. The analysis can only proceed if one chooses the sum of the ex- 
ponents to be constrained to 2/i + 1, whence the hypergeometric function 
becomes a simple algebraic function. Alternatively one can make a confor- 
mal transformation on the integral representing the graph of the left side of 
fig. 1. If z is the location of the integration vertex and the origin is chosen 
to be at the external end of the bosonic line then one can make the following 
change of variables, [fy|, 

*m — > % (3-1) 

which represents a conformal transformation. The part of the numerator 
involving fermion propagators transforms as, say, 

" ^fc^ (3-2) 

This results in a new integral of the same structure but where the exponents 
of the lines are related to the exponents of the original graph. In particular 
the exponent of the bosonic line becomes 2/i + 1 — J2i Pi- Therefore choosing 
this to be zero one obtains a simple chain integral and the result of fig. 
1 emerges naturally. We have detailed this type of conformal change of 
variables (3.1) and (3.2) here since it is fundamental to the problem and is 
the starting technique to the computation of all our graphs. 

With this basic rule one can define the conformal triangle construction 
for the vertex of the model. It is illustrated in fig. 2 for the cn/iip vertex of 
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(2.1) where the internal exponents a and b are defined to be such that each 
vertex of the triangle is unique ie 



2a + f3 = a + a + b = 2^ + 1 (3.3) 

It is worth noting that the way the model is defined 2a + (5 = 2/x + 1 — 2A 
so that the vertex of (2.1) is only unique at leading order in 1/N, since A = 
0(1/ N). The idea then of (6|, ||] is to replace all vertices, which are therefore 
not conformal, by a triangle which after either conformal transformations 
or applying the rule of fig. 1 yields the original vertex, up to a (amplitude) 
factor. In the conformal bootstrap approach all vertices of a Feynman dia- 
gram are replaced by conformal triangles. In other words each vertex of this 
new graph is unique or conformal and it is therefore elementary to observe 
that such graphs are then proportional to the basic (bare) graph by using 
conformal transformations which in turn allows one to write down simple 
equations for the Dyson equations of the Green's functions. 

For example, the 3-point vertex of (2.1) is given by an infinite sum of 
graphs, where we will use 1/N as our ordering parameter, and we have 
illustrated the first few graphs in its expansion in fig. 3 where the labels 
Ti will be used for later purposes. Each vertex of Tj and the left side are 
replaced by the construction of fig. 2 and the values of the integral they 
correspond to will be a function of the exponents a and (3 and the quantity 
z defined as z = fA 2 B where A and B are our earlier amplitudes and 
/ is the amplitude of the triangle. These are the three basic variables of 
the formalism and to deduce information on each we need three bootstrap 
equations. The vertex expansion of fig. 3 provides the first. Since each 
graph of fig. 3 will be a function of a, j3 and z and also by our arguments 
be proportional to the basic vertex, then if we denote by V the sum of all 
contributions from the graphs of fig. 3 it must be unity. Thus the first 
bootstrap equation is, || p~5f] , 

V(z,a,(3;0,0) = 1 (3.4) 

where we have displayed the arguments of the vertex function explicitly. 
The origin of the final two arguments will become apparent later when we 
have to introduce two regularizing parameters. Thus (3.4) represents the 
sum of all the contributions of the graphs in fig. 3. 



The remaining two bootstrap equations have been derived in [15| by 
following the analogous construction of || for a (f> 3 theory in arbitrary di- 
mensions. Rather than repeat that derivation we will record their form and 
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then make several comments. We have, [p|, 



d 



r(a - 1) = zt — V(z,a,(3;e,e' 
~N~ = -fe V \ z > a >P'> e > e ) 



e = e' = 



e = e'=0 



where 



TT^a 2 (a - l)a 2 {a - l)a{b)a((3) 



(3.5) 
(3.6) 

(3.7) 



T(p)(a-l) 2 (a-l) 2 a((3-b) 

The quantities e and e' are infinitesimal regularizing parameters introduced 
into the formalism by shifting a and (3 respectively by a — > a + 2e' , [3 — > (3 
+ 2e. Their introduction is necessary to avoid ill defined quantities in the 
derivation of the conformal bootstrap equations of the 2-point functions of 
(3.5) and (3.6). Basically a divergent integral arises which needs to be reg- 
ularized and this is achieved by the above shift though the integral appears 
multiplied by either e or e' so that only the residue of the simple pole, ie t, 
is required. 

The effect of the introduction of the regularization is that the sum of 
graphs contributing to the vertex function V become additionally functions 
of e and e', ie V(a, ft, z;e, e'). The graphs remain conformal because the 
internal exponents of the conformal triangle are adjusted (symmetrically) 
to preserve the uniqueness of each vertex. Thus with the regulators the 
vertex function is a sum of values which now depend on z, a, (3, e and e' 
and it is this function which appears in (3.4)-(3.6). Consequently, one has 
to compute the graphs of fig. 3 in the presence of the shift. 

With (3.5) and (3.6), it is now possible to deduce a master equation 
which determines 773. First, we let V\ denote the contribution to the vertex 
function from the one loop graph of fig. 3 and V2 the higher order pieces 
which will be 0(1/N) relative to V\. The former will later be expanded in 
powers of 1/N too. Then taking the quotient of (3.5) and (3.6) 



Nr(a - 1) 



dV_ 
de 7 



dV 



(3.8) 



On the left side of (3.8) we have a function of the exponents which can be 
expanded in powers of 1/N and will involve the anomalous exponents 77 and 
2A. Expanding to 0(1/N 2 ) the only unknown there is is 773 since rji, 7/2, Ai 
and A2 have already been determined. Thus to deduce 773 the right side of 
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(3.8) must be expanded to 0(1/N 2 ) and this will therefore involve the set 
of graphs contributing to V 2 . 

However, it is worth remarking on the functional structure of the graphs 
themselves. In representing each vertex by a conformal triangle one is re- 
placing a vertex of the theory which is non-unique by a construction where 
each vertex is unique. The non-uniqueness of this is reflected in the form 
of the value of the conformal triangle, which can easily be deduced by the 
method of subtractions discussed in [T^]. In particular it will be of the 
form 1/A where the residue is not important for the moment, [p9[| . In other 
words the variation from the value of the exponents which make the vertex 
conformal or unique appears as a pole in the value of the vertex. Therefore, 
a graph built out of n conformal triangles will have the structure h(A)/A n 
where h(A) is an analytic function of A. In the regularized theory the vari- 
ation from the conformal value is (A — e) and (A — e') in the vertices which 
are regularized. Therefore the form of a graph contributing to the vertex 
functions V x and V 2 will become h(A, e, e')/[A n_2 (A - e)(A - e')]. It is this 
structure which we now consider in the context of (3.8). Differentiating with 
respect to either regularizing parameters involves differentiating the residue 
or the simple pole in (A — e) or (A — e'). In the case of the latter the sum of 
all such contributions will be V/(A — e)\ or V/ (A — e')|, whilst in the former 
case it will be the sum of the residue contributions. The upshot of this is 
that within (3.8) when e and e' are set to zero the contribution from the 
poles cancels and only the residue needs to be considered. In other words, 

(3.9) 



Nr(a - 1) 



dV 



res. 



dV 

1 + A s7 



res. 



exactly where we have used (3.4) and the subscript res denotes differentiation 
with respect to e or e' of the residue of the vertex function. Therefore it is 
the expansion of each of these quantities to 0(1/N 2 ) which is required for 
773. More concretely we have 



Nr(a - 1) 

W) 



Ai 
1 + — 

N 



10 



de' 



10 



+ Ai 



A 1 



(dVi 
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V de' 



11 



de 



A 



1 

+ N2 


~A 2 ( dv ™ 

V de' 


dV 10 
de 


2 8V W 


(dV 10 




dV 10 




1 de 


V de' 




de 


) 



fdv 2 


dV 2 




\ de' 


de 


)] 



(3.10) 



where V\ = V\q + V\\/N and we now omit any comment on the evaluation 
symbol and z 2 is deduced from (3.4). Thus (3.10) forms the master equation 
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to deduce r/3. It has a two part structure. First, one requires the leading 
order graph V\ to be expanded to 0(1/N 2 ) and to be computed for non-zero e 
and e', whilst the graphs for V2 must be calculated. In |l5[ which established 
the bootstrap equations the vertex function for V\ was determined in order 
to check that the already known 0(1/N 2 ) results of ]l2| could be correctly 
determined. Moreover, the function was given at 0(1/N 2 ) for non-zero e 
and e' and we recall 



Q 3 



A(A-e)(A-e 



-exp[F(e,e',A)] 



(3.11) 



where 



F(e,e',A) 



+ 
+ 



5B/3 — 2B a -i 
Bq — -Bq-i — 

C a -1 



3Bn 



1 



a 



-\A-(B P -B )e 



a — 1 
7C/3 3Cq 



1 



Cp + Co — 2C a _i + 

Co — Cy3 — 2Cq,_1 + 



(a -l) 2 
2 



(a - l) 2 
2 



A 2 
Ae 
Ae' 



(3.12) 



with B x = tp(n - x) + 4>(x), Bq = -0(1) + ip(fJ>), C x = tp'(x) - - x), 
C = fj/(jj) - V'(l) and 



Q 



TT 2 ^a 2 (a-l)a((3) 

(«-i) 2 r(^) 



(3.13) 



We have illustrated the full T\ graph in terms of conformal triangles in fig. 
4 to show where the regulators appear in the appropriate exponents which 
is important for writing down the graphs of V2. It is crucial to note that in 
(3.12) one should substitute a = fi + \r] and (3 = 1 — 77 since the calculation 
of the graph was simplified by having the exponent A of the a line displayed 
explicitly. With (3.12) it is a trivial matter to deduce 772 agrees with (2.9) 
where the variable z is eliminated through 



A 3 



(3.14) 



at leading order from (3.4). 
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4 Computation of higher order graphs. 



All that remains is the evaluation of the higher order graphs which is far from 
a trivial exercise. The basic tools to compute the graphs are the conformal 
transformations (3.1) and (3.2) on each of the vertices of integration of the 
graph with the vertex joining to the external a line chosen to be the origin 



and the uniqueness rule of fig. 1 and its bosonic counterpart, (2^, 10]. Whilst 
it is possible to compute each graph with these methods easily when e = e' 
= which is required for the subsequent correction to (3.14), it turned out 
that it was not possible to evaluate the graphs exactly when e / 0, e' / 0. 
It transpired, however, that the difference, ie 



de' 



de 



(4.1) 



could be determined in each case, i = 2, ... ,5. There are four higher order 
graphs to consider. Whilst I^ and T§ are equivalent in the absence of the 
regulators, in the regularized version they are not equal since the bottom 
right external vertex of each Tj is regularized. In the evaluation of Tj we 
replace each unregularized conformal triangle with the original vertex and 
multiply the graph by —Q/A from (3.11), for each such vertex. It is only 
the regularized vertices which will give the contributions to the residue of 
V2 which are relevant for 7/3. 

Before discussing the determination of I^, and T§ we detail the cal- 
culation of T2 as it will illustrate some elementary steps which occur in all 
cases. Making a conformal transformation on T2 in the manner indicated 
yields a three loop integral which is illustrated in fig. 5 and we have evalu- 
ated two elementary chain integrals en route whose integration rule has been 
given in [p^ j, for example. However, this integral has two unique vertices 
which leads to the 2-loop integral of fig. 5 where we have now set A = 
since it does not contribute to the residue of T2 at this order in 1/N. (The 
poles in A emerge from factors such as a((i — A) which arise from a chain 
integral and the other from the integration of the two unique vertices.) In 



the notation of p7| it is (a, [3 — e' , a — e', /3, e + e'} but it cannot be evaluated 
by uniqueness methods for arbitrary e and e'. However, if we recall that it 
is the difference of the derivatives of the residues which is relevant for r/3 we 
can make use of the observation that 

■ d_ _ d_ 

. de' de 



d d 

— --)(a,P-e',a-e',P,e + e') 
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= ^7 (a,p-e',a-e',p,0) 



e'=0 



(4.2) 



to evaluate the relevant piece. The latter integral is equivalent to chain 
integrals and is simply 



v 2 (a -I, a -I, (3) 



Br, 



(a-1) 



(4.3) 



The overall value of the contribution is then given by carrying out the dif- 
ferentiation with respect to the regularizing parameters on the functions 
obtained after the successive integrations, which is elementary. We find 



(dT 2 


dT 2 \ 


\ de' 


de J 



2Q b v 2 {a-l,a-l,(3) 
(a - 1) 4 A 5 



Bp - B, 



1 



a-1 



(a-1) 



(4.4) 

Including the amplitude factor z 2 associated with the graph and using (3.14) 
we have 



(dV 2 , 2 dV 2 , 2 
V de' de 



(4.5) 



(/x- 1) 3 N 

in an obvious notation where we have substituted for the leading order values 
of a and (3 with respect to 1/N. 

The treatment of 1^ is equally as straightforward as that of F 2 and we 
note that 



(dT 3 


dT 3 \ 


\ de' 


de J 



7r 4 ^Q 7 a(2/i - 2)a(n - 1) 
A 7 (^-l) 6 

{2ji 2 -5^ + 4)* 



[(M - 1)(2a* - 1) 



5(2^ 2 -5^ + 4) 
207^1) 



(4.6) 



from which we obtain 

'dV 2 ,3 dV 2 ,3 



de' 



de 



(ji- 1) 2 N 



[( M -l)(2/x-l) 



(2ii 2 - 5n + 4)* 



5(2,u 2 -5/x + 4) 



2(/i-l) 

This calculation made use of the elementary chain integral 

- fM-fiU - j) [2(fi-a){n-P) + 3n-a-0\ 



(4.7) 



L 



(y 2 ) a ((x-y) 2 )P 



(a -1)09-1) 
x i/(a-l,/?-l,2/i-a-/0 + 2)(4.8) 
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The analysis of both and T5 turned out to be extremely intricate 
and involved and it is worth discussing several of the intermediate steps in 
one due to the presence of a common difficult integral in each. We consider 
T4 which after various transformations and elementary integrals leads us 
to consider the derivative with respect to 5 of the two loop integral of fig. 
6 where there is a trace over the lower right and one of the top fermion 
propagators. This is, in fact, a remnant of the fermion loop of the original 
graph. To evaluate this graph to 0(6) we restricted attention to the leading 
order values since that leads to an integral where most of the exponents are 
unity which is convenient for taking a fermion trace. For instance, taking 
the trace over the open fermion propagators and dividing by the appropriate 
factor trl, fig. 6 is equivalent to 

±tr[(l, 2, 1, 1 - 5, m - 1) - (0, 2, 1, 1 - 6, fx - 1) - (1, 1, 1, 1 - 5, fi - 1)] (4.9) 

The second term is trivial to deduce whilst the third becomes trivial after 
making the transformation <— in the notation of flQ|| . The hardest part of 
the integral lurks within the first term, which after performing the trace 
explicitly yields one trivial graph and 



(1,2,1,1 -5,ii-2) - (1,1,1,1 -S it i-1) 



(4.10) 



The 5 expansion of an integral similar to the second term of (4.10) arises 
in the analogous 0(1 /N 3 ) calculation in the pioneering work of ||. Unfor- 
tunately it is not possible to give a closed form for the 0(5) correction in 
terms of elementary functions like \& and for all dimensions. Instead one 
is forced to leave the 0(5) correction defined as the quantity I(fi). Thus 

1-5,^-1, n-l + S) = ChT(l,l)[l + 51(h) + 0(5 2 )} (4.11) 

where we have used the result ChT(l, 1 — 5) = ChT(l, 1+5) and the function 
ChT(a,(3) is defined to be (a, fj, — 1, /J, — 1, (3, fj, — 1) for all a and (3 in the 



notation of [ 10 1 where it was evaluated exactly as 



ChT(a,(3) 



vr 2 ^a(2// - 2) 

rV-i) 



a(a)a(2 



+ 



a(P)a(2 - (3) 



(1 - (3)(a + (3 - 2) (1 - a)(a + (3 - 2) 



+ 



a(a + (3 — l)a(3 — a — (3) 



(a-l)(P-l) 



(4-12) 



Whilst a closed form for I(fi) is not available it can be analysed using the 
Gegenbauer polynomial techniques of p0| to obtain a set of double sums 
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over T-functions and its e-expansion near d = 2 + e can be given to several 
orders. Then with the definition 



we have 



3(m- 
C(3)6 2 C(4)6 3 



+ 



■ + m 

13C(5)e 4 



+ 0(e 5 



(4.13) 



(4.14) 



6 8 48 

where the first few terms of (4.14) were given in [10| and later ones in pl| , [32] ], 
It turns out that in three dimensions an exact expression can be deduced as 



2 In 2 + 



2vr 2 



(4.15) 



The remaining integral of (4.10) also contains which can be deduced 
by various transformations of bosonic two loop integrals given in [PJ and 
recursion relations of One result which was required in this and which 
is worth recording is 



— <M-l + e,2,l,M-l,At-2> 



+ 



7r^a(l)a(2// - 2) 

(2/j ~ 3) 
2( M -l)(/i-2) 



2) 2 
1 



+ 



(2/x-3)(//-3) 2(^-1) (2/*-3)* 



(m-i) 



( M - 2) (/x - 2) 
- 3(// - 2)I(p) (e + 



1 



01-1)2 

(4.16) 



Consequently we have that the coefficient of the 0(<5) term of (4.9) is 



vr 2 ^a(l)a(2^ - 2) 



2(M - 1) 
+ Qx-5)0 + 



3(/x- 1)03 + 



(m-i) 



(2// - 3) 



+ 



(m-i) (^-i) 2 



(4.17) 



Collecting all the contributions to the integral and evaluating the derivative 
of the residue gives the relatively simple result 



2.1 



2.1 



8e' 



N 



303 



+ 



33 



+ 



2(^-1)2 



( M -8)0 (2/X-3) 
2( M -1) + (^-l) 3 



(4.18) 



16 



The procedure for T$ is equally as tedious but does not merit extensive 
coverage since the key techniques have already been covered in the discussion 
of IV A similar integral to fig. 6 occurs but this time the regulator is not 
on the line with exponent (5 but on the other bosonic line appearing as 
(// — /? — e). Its O(e) contribution is 



n 2 »a(l)a(2fi - 2) 

2(0^1) 



3(0-1)6- 



3~ 



(M-l) 



(3/i + 2)e + 



2^ 



(M-l) 



(4.19) 



as an intermediate check on the calculation. Overall we find 



2,5 



2,5 



de 



N 



36H 



+ 



3~ 



+ 



0* 



46 



2(0-1)2 (/x _l)2 

(2 M -5)(Ai-2)" 



(M-l) 2(0-1)3 



(4.20) 



This completes the discussion of the contributions to V2. However, it is 
worth recording the values of each of the Tj when e = e' = since they are 
required for the higher order corrections to (3.14) which determines Z2- We 
found that 

Q 5 i/ 2 (a-l,a-l,/3) 



A 5 



(a-1) 



(4.21) 



r 3 = 
r 4 



Q 7 u s (a -l,a-l,0) 



A 7 

r 5 



i/(/J,/?,2 A *-2/J)[2( M -/J)( M -/J-l)+ A *] (4.22) 

(2/X-3)" 



(a - l) 6 /3 2 

7r 2 ^a(l)a(2/i - 2)QV(a - 1, a - 1, /3) 



2A 7 (a - 1) 



36 



I) 2 . 
(4.23) 



5 Discussion. 

With the explicit values of the V2 contributions to (3.10) calculated it is now 
possible to deduce 773 although the remainder of the calculation involves a 
substantial amount of tedious algebra. We find 



+ 



06 



(3// - 1) , 0(0 - 16) 



(20 - I] 



4(0 
2$ 



+ 



I) 2 4(0 
3/i 2 H 



1) 



1 



6// 



(m-i) ; 



+ 



30 2 6~ 



^ 3/i 2 9^ 3(2// 



1)2*2 



2( M -l) 2 2(0 - 1)3 2(»-l) (0-1) 2(0-1 
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, T , ( M 3 M2/X-1) ^(2^-5^ + 4) _ 3 
2(/x-l)3 + 2(^-1)3 ^2/x 
35 18 9 \ 1 3 7 

+ 7 7777 + 777 TTq I +^"4 



2(m-1) (/j-1) 2 2( M -l)3y 2 M 2 M ( M _i) 

(5.1) 



+ 13 | 8 | 2 ^(2/i- l)(/x-3) 



2( M -1)2 ( M -l)3 ( m -1)4 ( M _l)3 

which is the main result of this paper and represents the first 0(1/N 3 ) 
analysis in this model in arbitrary dimensions. (Previous 0(l/iV 3 ) analysis 
was strictly two dimensional and examined the corrections to the a field 
mass, pi].) We can evaluate (5.1) in three dimensions to obtain 



256 

713 = 27^ 



47vr 2 2l 189C(3 167 

+ 9vr 2 n 2 

12 2 9 



(5.2) 



which has a similar structure to the analogous quantity in the O(N) bosonic 
a model in terms of the appearance of numbers such as In 2 and C(3). 

We conclude with several remarks. One important point of the analy- 
sis we have described here is the adaptation of conformal techniques and 
the conformal bootstrap programme to a model with fermions. Whilst the 
interaction is relatively simple it ought now to be possible to perform an 
analogous calculation in a gauge theory such as QED where the only essen- 
tial difference is the appearance of a 7-matrix at each vertex. This is partly 
motivated by the fact that 772 could be deduced via the self consistency 
approach of |L0| in [g^] and it ought therefore to be possible to construct 
analogous bootstrap equations for that model to allow it to be probed as far 
as is analytically possible. 
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Note added. Whilst in the final part of this work we received a preprint, 
p6|, where 773 is also recorded and we note that both results are in agreement. 
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Figure Captions. 

Fig. 1. Uniqueness rule for a fermionic vertex. 

Fig. 2. Conformal triangle. 

Fig. 3. Expansion of 3- vertex. 

Fig. 4. Regularized 3-vertex graph. 

Fig. 5. Intermediate graphs in the evaluation of T2. 

Fig. 6. T4 after several integrations. 
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